Neutrino mass hierarchy extraction using atmospheric neutrinos in ice 
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We show that the measurements of 10 GeV atmospheric neutrinos by an upcoming array of 
densely-packed phototubes buried deep inside the IceCube detector at the South Pole can be used 
to determine the neutrino mass hierarchy for values of sin 2 2(9i3 close to the present bound, if the 
hierarchy is normal. These results are obtained for an exposure of 100 Mton years and systematic 
uncertainties up to 10%. 
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I. INTRODUCTION 

Neutrino physics has undergone a true revolution over 
the last decade, with a large number of experiments of 
different types providing evidence for neutrino oscilla- 
tions and thus for physics beyond the Standard Model. 
This has led to new insights into the possibilities for 
new physics, but while some old questions have been an- 
swered, many new questions have emerged 

An important question that still needs to be resolved is 
the determination of the neutrino mass hierarchy. A very 
large effort has been dedicated to planning new experi- 
ments for measuring all neutrino oscillation parameters 
2]. The mass hierarchy can be determined using mat- 
ter effects on oscillations inside the Earth. This however 
requires a long baseline, a very large detector and an in- 
tense beam. In addition, parameter degeneracies have to 
be resolved using a combination of experiments. 

Cosmic ray interactions in the atmosphere give a nat- 
ural beam of neutrinos. These atmospheric neutrinos in 
the GeV range have been used by the Super-Kamiokandc 
detector to provide evidence for neutrino oscillations. 
The large size of neutrino telescopes such as AMANDA, 
IceCube and KM3NeT makes possible the detection of 
a large number of atmospheric neutrino events with a 
higher energy threshold, ~ 100 GeV, even though the 
neutrino flux decreases rapidly with energy (~ E~ 3 ). 
Built to detect neutrinos from astrophysical sources, dark 
matter annihilation, etc. these ice/water Cherenkov 
detectors typically have high detection threshold energy 
to avoid the large background from atmospheric neutri- 
nos. Since neutrino oscillation effects are quite small at 
high energies, high energy atmospheric neutrinos provide 
little information about important issues such as neutrino 
mass hierarchy, mixing angles, etc. 

Recently, a low energy extension of the IceCube de- 
tector has been planned [4J. Its goal is to significantly 
improve the atmospheric muon rejection and to extend 
the IceCube neutrino detection capabilities in the low en- 
ergy domain, possibly to muon energies as low as 5 GeV, 
depending on the density of the photo tubes. The pro- 
posed instrumented volume is 5,000-10,000 kton by con- 



servative estimates. Such a low threshold array buried 
deep inside IceCube will open up a new energy window 
on the universe. It will search for neutrinos from sources 
in the Southern hemisphere, in particular from the galac- 
tic center region, as well as for neutrinos from WIMP an- 
nihilation, as originally motivated. Here we provide an 
additional and independent motivation for building such 
an array, namely to explore neutrino oscillation physics. 

In this article we show that the deep core extension of 
IceCube provides a great opportunity for detailed oscilla- 
tions studies of atmospheric neutrinos and makes possible 
the determination of the neutrino mass hierarchy. This 
is extremely important given that long baseline experi- 
ments with comparable sensitivity might take a very long 
time to build and collect data, while the IceCube deep 
core will be built in the near future and will accumulate 
high statistics relatively fast. Our study indicates that 
for a total mass of the instrumented volume times ex- 
posure of 100 Mt yr (roughly equivalent to 20 years of 
running a 5,000 kton detector, or 10 years of running a 
10,000 kton detector in an optimistic scenario), neutrino 
mass hierarchy can be determined at least with 90% con- 
fidence level assuming the current best-fit values of the 
oscillation parameters, and for values of #13 close to the 
present bound. 

We start in section |TT] with a brief review of what we 
know about neutrino oscillation parameters and what we 
expect to learn in the near future. We also discuss matter 
effects inside the Earth, which play an important role in 
the analysis. In section ITTT1 we describe our main analysis 
and in section lTVl we discuss the backgrounds and system- 
atic uncertainties that affect this analysis. We present 
the main results in section IVl and discuss them in section 
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II. NEUTRINO OSCILLATIONS 

Neutrino data from solar, atmospheric, reactor and 
accelerator experiments is well understood in terms of 
three-flavour neutrino oscillations. Two Am 2 values and 
two (large) mixing angles are well determined, while the 
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third mixing angle is limited to be very small. The CP- 
violating phase (5) is completely unconstrained. In ad- 
dition, the sign of Am^ is also unknown and will be 
the focus of our study. The two possibilities, Am^ > 
or Amjj-L < correspond to two types of neutrino mass 
ordering, normal hierarchy and inverted hierarchy. 

The best fit oscillation parameter values obtained from 
present data are @, Q : 
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sin 2 26* 23 
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with 99% CL allowed regions given by: 
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and sin 2 26»i 3 < 0.15 for Am§ x = 2.5 x 10" 3 eV 2 . No- 
tice that an extra unknown in the neutrino oscillation 
scenario is the octant in which 6*23 lies, if sin 2 26*23 7^ 1- 
This has been dubbed in the literature as the 6*23 octant 
ambiguity. 

In the near future, long baseline experiments like MI- 
NOS and OPERA will improve the current precision on 
Am 2 ^ and possibly discover a non-zero value of #13, if this 
is close to the present upper limit. In a few years, reactor 
experiments like DoubleChooz and DayaBay will provide 
improved sensitivity to 6*13 . This information can be used 
as input in our analysis, reducing some of the parameter 
uncertainties. 

In the past, atmospheric neutrinos in the Super- 
Kamiokande detector have provided evidence for neu- 
trino oscillations and the first measurements of lAmJjJ 
and sin 2 26*23. It is also known that, covering a large 
range of energies and pathlengths and using matter ef- 
fects inside the Earth, they can be in principle sensitive 
to sub-dominant neutrino oscillation effects like 6*13 and 
the mass hierarchy 0J^, i, ©, M, [U E El • 



Matter effects [19J, [20j in long baseline and atmospheric 
neutrino oscillation experiments depend on the size of the 
mixing angle 6*13 which governs the transitions v e <-> v^. T 
driven by the atmospheric mass squared difference A31 = 
Am| 1 /2£J. The effective 6*13 mixing angle in matter in a 
two-flavour framework is given by: 



sin 2 26*™ = 



sin 2 26 



13 



sin 2 26»i3 



cos 26*13 T 



s/2G F N e 
A31 



(9) 



where the minus (plus) sign refers to neutrinos (antincu- 
trinos), N e is the electron number density in the medium, 
V2G F N e (eV)= 7.6 x l()- 14 Y e p (g/cm 3 ) and Y e , p the 
electron fraction and the density of the medium, the 
Earth interior in our case. The electron fraction Y e is 



0.466 (0.494) in the core (mantle) and we follow the 
PREM [2l[ model for the Earth's density profile. Equa- 
tion © implies that, in the presence of matter effects, 
the neutrino (antineutrino) oscillation probability gets 
enhanced if the hierarchy is normal (inverted). Making 
use of the different matter effects for neutrinos and an- 
tineutrinos seems therefore the ideal way to distinguish 
among the two possibilities: normal versus inverted mass 
hierarchy. Matter effects are expected to be important 
when the resonance condition: 



Affij! cos (26*13) 



2V2G F N e E 



(10) 



is satisfied. The precise location of the resonance will de- 
pend on both the neutrino path and the neutrino energy. 



For Ar 



'31 



2.5 x 10 3 eV 2 and distances of several thou- 



sand kilometers the resonance effect is expected to take 
place for neutrino energies O(10) GeV. The pathlength 
traveled by atmospheric neutrinos is: 



L(c„) = R®(^(l + l/R®y 



(11) 



where R^ — 6371 km is the radius of the Earth and 
I ~ 15 km is the typical height at which neutrinos get 
produced in the atmosphere. The cosine and sine of the 
nadir angle of the incident neutrino are denoted by c„ 
and s„, respectively. Since upward going (c„ — * —1) at- 
mospheric neutrinos traverse the dense core of the Earth, 
they provide an excellent tool to tackle the neutrino mass 
ordering. 

Indeed the idea of exploiting matter effects in atmo- 
spheric neutrino oscillations to distinguish the type of 
hierarchy has been extensively explored in the litera- 
ture d, @, H, [l(| ■ In general, the former studies exploit 
muon calorimeter detectors, such as MONOLITH [l6| . 
MINOS [l7| or INO [lj] in which the muon charge can 
be determined. The measurement of the number of pos- 
itive and negative muons in the 1 — 10 GeV energy re- 
gion allows then for a direct extraction of the neutrino 
mass hierarchy, simply by looking in which channel (neu- 
trino or antineutrino) the signal, via matter effects, is en- 
hanced. However, it has been pointed out, and carefully 
explored, that the detection of atmospheric neutrinos 
which have crossed the Earth by future planned mega- 
ton water Cherenkov detectors could also determine the 
neutrino mass hierarchy, provided the mixing parameter 
sin 2 26*13 is not very small (Til . [T2I [l3j , even when these 
detectors do not allow for a charge discrimination of the 
leptons. The detection of low energy neutrinos in a higher 
density photo-multiplier array within the IceCube instru- 
mented detector volume opens up the possibility of not 
only exploring the atmospheric oscillation pattern [l5l |. 
but also exploring how well the neutrino mass hierarchy 
could be measured in the largest water/ice Cherenkov 
detector available in the near future. In the next section 
we present the details of the analysis proposed here. For 
our numerical analysis, unless otherwise stated, we will 
use the best fit values quoted earlier in this section for 
the oscillation parameters. 
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III. ANALYSIS 

The IceCube detector has the ability to measure sep- 
arately the muon tracks and electron/tau generated cas- 
cades, thus providing good flavour identification in some 
energy ranges. This would be extremely useful for an 
oscillation analysis, especially when searching for sub- 
dominant effects. In particular, the atmospheric neutrino 
sensitivity to the mass hierarchy comes from the matter 



effects on 



oscillations and previous 



studies dealing with water Cherenkov detectors have used 



J 



the electron signal to extract this information. In the low 
energy range relevant for neutrino oscillations however, 
it is extremely hard in IceCube (and even in the deep 
core array) to obtain information about neutrino direc- 
tion and energy for electron cascades. 

We thus focus here on the /x-like contained events pro- 
duced by the interactions of atmospheric upward going 
neutrinos in deep ice. Formally, the expected number 
of muon neutrino-induced contained events in the i- and 
j-th energy and cosine of the nadir angle (c„) bins read: 



2nN T t 

Vdct 



Ei+Ai 



dE v 



dc v V,, 



' 1 ^ ^K)^>e)^ + 



'"(.("=) cc 



dE v dVl 



r 



(12) 



where A* and Aj are respectively the energy and c„ bin 
widths, is the number of available targets, Vdet is the 
total volume of the detector, t is the exposure time, d<^„'s 
are the atmospheric (anti)neutrino differential spectra, 
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cr cc is the CC (anti)neutrino cross section and is the 
effective detector volume. For a detector with cylindrical 
shape of radius r and height h, is given by [151 ] 



V li (E^6)=2hr'zxcsm\ VI- 




(13) 



r 



where R Pl (E fJi ) is the energy-dependent muon range in ice. 
For a cc we use the charged current (anti)neutrino inter- 
action cross-sections in 25]. It is useful to note that in 
the relevant energy range the anti-neutrino cross-section 
is smaller than the neutrino cross-section by about a fac- 
tor of two. This difference is what allows for the (statisti- 
cal) discrimination between neutrinos and anti-neutrinos 
and thus between normal and inverted hierarchy in this 
experiment. 

Equation Q12[) contains the atmospheric electron and 
muon (anti) neutrino fluxes, f^"^Q ■ For the results pre- 
sented in this study we use the results from Refs. [22] ]. 
The atmospheric neutrino fluxes from Refs. [23] have 
also been used, and, overall, we obtain a similar dif- 
ference between the number of muon neutrino induced 
events for the normal and inverted hierarchies both in 
energy and c„. The absolute electron and muon atmo- 
spheric (anti)neutrino fluxes are found to have errors of 
10% — 15% in the energy region of interest here [24j |. 
Those errors are mostly induced by our ignorance in 
modeling hadron production, although the situation is 
expected to improve with HARP and MIPP data. The 
uncertainties quoted above are reduced for the neutrino- 
antineutrino flavor ratio case [30| . i.e. for v^jv^ and 



v e /v ei where the uncertainty is ~ 7% in the energy range 
we explore in the current study. Even smaller uncertain- 
ties are expected when the muon-to-electron flavor ratio 
(z/p + 9^)/{v e + v e ) is considered. We will comment on 
the impact of the atmospheric neutrino flux uncertainties 
on our results below, when including systematic uncer- 
tainties to our numerical analysis. 

The energy of secondary muons from CC interaction 
in the 10-100 GeV neutrino energy range of interest here 
is (E^) — 0.52 E v and 0.66 E v , respectively for neutri- 
nos and antineutrinos [2||. We illustrate in Fig. [1] the 
expected /i-like contained events in 5 GeV muon energy 
bins for a combined detector mass times exposure of 50 
Mt yr. From left to right, the panels depict the con- 
tained /i-like events within the (—1,-0.9), (—0.9,-0.8) 
and (—0.8, —0.7) c v bins. We assume sin 2 2$i 3 = 0.1 and 
S = along with other best-fit parameters in Eq. ©. 
Although we used a detector geometry of 1 km height 
and ~ 40 m radius, all events are contained in these c v 
bins except for the highest energy bin: (25,30) GeV. As 
we will shortly see, the oscillation signals solely affect the 
low energy events. Thus our results are valid for a va- 
riety of detector geometries, only affected by the total 
instrumented volume times the exposure. 

For the first c v bin used in Fig. [T] the (anti)neutrinos 
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FIG. 1: From left to right, number of contained /i-like events in a 5000 kton detector after 10 years exposure, within the 
( — 1,-0.9), (—0.9,-0.8) and (—0.8,-0.7) c v bin, assuming sin 2 2#i3 = 0.1, 5 = and an energy bin size of 5 GeV, with a 
muon energy threshold for detection of 5 GeV. We show the la statistical errors. The blue crosses (red circles) denote normal 
(inverted) hierarchy. 



are almost vertically upward going and in their way to- 
wards the detector they have crossed the high density 
Earth core. For the other bins (anti)neutrinos still tra- 
verse a significant amount of matter through the Earth. 
Note that a finer angular bin than Ac„ ~ 0.1 is not pos- 
sible because the reconstruction of primary neutrino di- 
rection is expected to be poor in this energy range due 
to the intrinsic spread in charged lepton-neutrino scat- 
tering angle. The resonance is expected to be located at 
low energies, see Eq. (|10[) . and the maximum difference 
between normal and inverted hierarchies is observed in 
the (5, 10) GeV energy bin. In higher energy bins, the 
effect is totally negligible. 

Figure [1] helps understanding the energy and angular 
bins that should be considered for the numerical anal- 
ysis developed in the current study. We will exploit 
exclusively the first three muon energy bins (i.e. 
within the (5,10) GeV, (10,15) GeV and (15,20) GeV 
energy ranges), and three angular bins, i.e., c„ within 
the (-1,-0.9), (-0.9,-0.8) and (-0.8,-0.7) ranges). It 
is precisely in those bins where the sensitivity to the neu- 
trino mass ordering is significant, and all the muon events 
are fully contained. Since the atmospheric neutrino spec- 
tra are very steep (~ E~ 3 in the relevant energy range) 
and, as we will discuss next, this energy range also corre- 
sponds to a neutrino oscillation maximum, most events 
are concentrated in the lowest energy bin. It is thus im- 
portant to keep the muon energy threshold for detection 
as low as possible (between 5 and 7 GeV) to collect large 
statistics as well as for extracting the neutrino mass or- 
dering. Note that, as discussed above, 5-7 GeV muon en- 
ergy threshold corresponds to approximately 10-15 GeV 
initial neutrino energy. With a proposed 7 m spacing 
between phototubes on the deep core strings, such a low 
energy threshold is achievable. We next discuss the os- 
cillation effects in the event rates in Fig. [TJ Notice that 
for the three angular bins illustrated here, the number of 



/i-like events in the first energy bin is larger for the in- 
verted hierarchy case than for the normal hierarchy case. 
The reason for that is due to the presence of matter ef- 
fects at these zenith angles for f M 's, which are the ones 
dominating the statistics. This can be well seen in Fig. [21 
which shows P„ e ^„ M (= P eM ) and P v ^ Vfl _ (= P^) for both 
normal and inverted hierarchy, for sin 2 20i3 = 0.1 and 
sin 2 29 13 = 0.06. 

If the hierarchy is normal (inverted), the P MM survival 
probabilities, in the ^10-20 GeV energy range for neutri- 
nos (which are responsible to produce 5-10 GeV muons), 
are suppressed (enhanced), due to matter effects, and 
therefore a smaller (larger) number of CC interac- 
tions are expected in the detector. If the matter density 
is constant, and the contribution from the solar terms is 
negligible, the P MAt survival probability is given by (see 
Refs. [H,[ll) 



p — 
1 mi. — 



1 - cos 2 0g sin 2 20 



1.27 



2:t 



(14) 



Am 2 ! + A + {Aml 1 ) m \ L 

E 



sin 2 6»™ 3 sin 2 2(9 23 



x sm 



1.27 



Am| x + A- {Aml 1 ) m \ L 



E 



sin 4 6>23 sin 2 29^ 3 sin 2 



1.27(Am 2 1 ) m | 



where A = 2V2G F N e E, 9f z is given by Eq. © and 
(AmL) m = Am 2 x 



sin 2 20i3 



cos 20i3 T 



2V2G F N e E 



Am 2 ! 



where the minus (plus) sign applies to neutrino (antineu- 
trino) flavor transitions. Due to the presence of resonant 
matter effects in the first angular and energy bins, P MM 
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can be very different for normal and inverted hierarchies 
(mostly due to changes in the first term in Eq. (fTS"|) ). 

As seen in Fig.[2j the v e — > probability is quite small 
in the relevant energy range (negligible for inverted hier- 
archy) . Its contribution to the final muon-like event rate 
is made even smaller by the fact that the atmospheric 
v e flux at these energies is much smaller than the 
flux. The experiment we are considering is thus mostly 
exploiting matter effects in the disappearance — > 
channel and therefore is in many ways complementary 
to the appearance experiments. While the matter effects 
are a small correction in the survival probability, they 
are sufficient to provide a difference between the differ- 
ent mass orderings because of the very large number of 
events. 

Note that in Fig. Q] the difference between event rates 
for the two hierarchies increases (although the overall 
rates decreases) for c v bins (—0.9, —0.8) and (—0.8, —0.7) 
compared to the (—1, —0.9) bin. This is because the res- 
onant matter density for neutrino energies in the first 
energy bin < E v >= 15 GcV is ~5 g/cm 3 which is lower 
than the densities that the neutrino crosses if c„ is in 
the (—1, —0.9) region, but gets closer to the ones in the 
shallower c„ region. 



IV. BACKGROUNDS AND SYSTEMATIC 
UNCERTAINTIES 

The main backgrounds to the signal we are exploiting 
in the current study are atmospheric downward going 
muons from the interactions of cosmic rays in the atmo- 
sphere and tau (anti)neutrinos from v^^i^^.e) — ► v T (v T ) 
transitions. The cosmic muon background can be elimi- 
nated by angular cuts and in the Ice Cube deep core is 
significantly reduced compared to the IceCube detector. 

The tau neutrino background can be included in the 
analysis as an additional source of fi-\ike events. Tau 
(anti)neutrinos resulting from atmospheric neutrino fla- 
vor transitions will produce a r lepton by CC interac- 
tions in the detector effective volume. The tau leptons 



produced have an ~ 18% probability of decaying through 
the t~ — > [i^v^Vr channel. 

The secondary muons can mimic muons from CC 
interactions and must be included in the oscillated signal. 
The energy of a v T needs to be about 2.5 times higher 
than a to produce, via tau decay, a muon of the same 
energy. But the atmospheric neutrino flux has a steeply 
falling spectrum, so one would expect this tau-induced 
muon background not to be very large. It is however sig- 
nificant (~ 10%) due to the fact that, as seen in Figure^ 
the first maximum in the — > v r oscillation probability 
(minimum in the — > survival probability) falls ex- 
actly in the energy range of interest and for a large range 
the v r flux can be significantly larger than the flux. 
These events significantly change the energy spectrum of 
the measured muon-like events and contain information 
about the main oscillation parameters, A 31 and # 2 3 • 




E v [GeVJ 

FIG. 3: survival probability and —> v T oscillation prob- 
ability for c„ = — 1, sin 2 28n = 0.1 

The uncertainties in the atmospheric neutrino flux 
have been discussed in the previous section and they af- 
fect the analysis. It is however possible to use the data 
itself to improve some of the errors introduced by these 
effects, by considering energy and angular bins where os- 
cillation effects are not important as a reference and thus 
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canceling out some of these uncertainties in the analysis 
(see also HI). 

The uncertainties in other oscillation parameters also 
affect the possibility of determining the neutrino mass hi- 
erarchy. We employ a full three-flavour oscillation analy- 
sis. It is however easy to understand that the solar Am 2 
and mixing angle have almost no contribution. This is 
due to the high energy values we are considering, for 
which even the largest distances traveled by atmospheric 
neutrinos do not allow the Am^iL/E term to become 
significant. The present uncertainties in lAmjj-J and #23 
(see Eq. (jHJ) will be improved by present accelerator ex- 
periments. In addition, the atmospheric neutrino data in 
the IceCube deep core can be used to extract these pa- 
rameters independently from the sub-dominant effects, 
by comparing different angular bins. The values of #13 
and 8 are the most uncertain and usually hardest to dis- 
entangle from the mass hierarchy. A measurement of 
the 0i3 mixing angle by the near future reactor experi- 
ments would guarantee the possibility of extracting the 
mass hierarchy from the atmospheric neutrino data in 
IceCube. We consider in our analysis a range of values 
for 6>i3 and 8 in order to assess the parameter space for 
which the neutrino mass hierarchy can be determined. It 
is important to emphasize however that the sensitivity 
to the CP-violating phase <5 and the value of 623 in this 
experiment is quite small due to the fact that the signal 
is dominated by the z/ M disappearance channel. 



V. RESULTS 

In this section we present our results in terms of a min- 
imum x 2 analysis in the (sin 2 2^13; 8) parameter space. 
For a particular hierarchy (h) and (sin 2 29i 3 ] 8) parame- 
ters chosen by nature, we consider the number of /i-like 
events N£* h (sin 2 2^13; 8) measured by an experiment in 
the i- and j-th muon energy and c v bins (see Eq. (fT2|) V 
These events include the Vu and P„ signal, as well as the 
background secondary muons from v T and Cv's. 

The x 2 statistics, for a "theoretical" model of hierarchy 
(hf) and parameters (sin 2 29' 13 ; 8') is then defined as 

X 2 ,(sin 2 20' 13 ;<n= £ £ 

i=l,3 j=l,3 

" iVg^sin 2 2fl 13 ; 8) - iV^ v (sin 2 20[ 3 ; 8') Y 
a- i (sin 2 20 13 ;<5) J 

Here N^ h ,(sin 2 29' 13 ; 8') is the expected event number 
from both signal i^'s and background ^ T 's given a "the- 
oretical" model. The variance cr™^ (sin 2 28i 3 ; 8) is cal- 
culated from experimental events with or without sys- 
tematic uncertainties. We minimize x 2 m Eq. <| 1 6j) for 
(sin 2 28' 13 ; 8') parameters (i.e. 2 d.o.f). When nature's 
choice or "true" hierarchy is h (normal for example), then 
the "wrong" theoretical model of hierarchy h' ^ h (in- 



verted in this case) is rejected if 

min (xl'Tth) - min (Xh<=h) > a ( 17 ) 

in the (sin 2 26»i 3 ; 8) parameter space. The 68%, 90%, 95% 
and 99% confidence levels (CL) are defined for a = 2.3, 
4.61, 5.99 and 9.21 respectively, for 2 d.o.f statistics [271 ]. 
Note that our choice of 2 d.o.f statistics is rather con- 
servative as explained in the Appendix of Ref. [Hj . Our 
90% CL translate into the 97% CL for 1 d.o.f statistics. 

Fig. [5J left (right) panel shows the % 2 results in the 
(sin 2 20i3 > 8) plane for a measurement of the hierarchy at 
the different confidence levels quoted above, exploiting 
the muon- like contained events in a 5, 000 kton detector 
after 20 years exposure (100 Mt yr), for normal (inverted) 
hierarchies. The sensitivity to the mass hierarchy is bet- 
ter in case nature has chosen the normal hierarchy than 
in the case for the inverted hierarchy. The reason for 
that is that in the normal mass ordering scenario, matter 
effects decrease (increase) the muon neutrino (antineu- 
trino) survival probability. In the inverted mass order- 
ing scenario, matter effects decrease (increase) the muon 
antincutrino (neutrino) survival probability. Notice the 
observable we are exploiting here contains both the at- 
mospheric neutrino and antineutrino fluxes, see Eq. \12\ . 
Since in the energy range of interest the neutrino fluxes 
are roughly twice the antineutrino ones, the overall ef- 
fect for the normal hierarchy is larger, i.e. no significant 
cancellation among the effects for neutrinos and antineu- 
trinos is present. In the inverted hierarchy case, there 
exists a partial cancellation and therefore the sensitivity 
results are not as promising as the ones for the normal 
hierarchy. 

The dependence of the x 2 results on the value of the 
CP-phase 8 is extremely mild, as expected, due to the 
fact that the muon neutrino survival probability in mat- 
ter Pfj.fi, Eq. (fl~5|) . does not depend on 8 in the limit of 
negligible solar effects. Even when solar effects are not 
negligible, P MM in matter depends exclusively on cos 8. 
Consequently, one would expect the same results for 8 
and —5, as can be clearly noticed from Figs. 2] Also, no- 
tice how the maximum (minimum) sensitivity is reached 
at 8 = 0° (8 = 180°), due to the change of sign in the 
cos 8 function. 

Figs. [5] show the equivalent to Figs. Q] but including 
in the x 2 analysis performed here a 10% overall system- 
atic error, which includes uncertainties from atmospheric 
neutrino fluxes and cross sections, among others. The im- 
pact of the systematic errors on the sensitivity/exclusion 
curves is significant, therefore, keeping systematic uncer- 
tainties below the 10% level is crucial for extracting the 
neutrino mass hierarchy. While this value might seem op- 
timistic given our present knowledge of fluxes and cross- 
sections, it is expected these uncertainties will be signifi- 
cantly smaller on the time scale relevant for the IceCube 
deep core experiment. This is likely to happen in two 
ways: a number of precise cross-section measurements in 
the relevant energy range will be performed in the near 
future, and, most importantly, the systematic uncertain- 
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FIG. 4: Rejection regions of the "wrong" hierarchy model in the (sin #13, 5 cp ) plane when the "true" hierarchy is normal (left 
panel) or inverted (right panel) as indicated in the heading of each plot. Different lines correspond to rejection regions of the 
"wrong" hierarchy at 68%, 90%, 95% and 99% CL (2 d.o.f.) using the muon-like contained events in a detector of mass times 
exposure of 100 Mt yr (5, 000 kton detector after 20 years of data taking, or 10, 000 kton detector after 10 years of data taking 
in the optimistic scaenario). We used the best fit parameter values in Eqs. (l)-(4). 
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FIG. 5: Same as Fig. [4] but including a 10% systematic error in our x 2 analysis. 



ties can improve with data from the Deep Core itself. 
Our analysis only makes use of the lowest energy range 
and almost straight up-going neutrinos, which are most 
sensitive to the signal we are interested in. Neutrinos 
from all other directions and at higher energies are no 
longer sensitive to the sub-dominant neutrino oscillation 
effects, but they contain important information about 
fluxes, cross-sections and potentially main oscillation pa- 
rameters. They can thus be used as a "reference" in order 
to minimize the systematic uncertainties for our analysis. 



We also explore the impact of the #23-octant ambigu- 
ity. Following the current 95% CL allowed limits for the 
mixing parameter sin 2 2#23, we illustrate here the results 
from our \ 2 analysis for #23 = 40° and #23 = 50°. Figs. [6] 
show the x 2 analysis results, including a 10% overall sys- 
tematic error for 823 = 40°. Figs. [7] show the equivalent 
but for 6» 23 = 50°. 

If #23 7^ 45°, the results for the normal hierarchy be- 
come slightly worse, since sin 2 2# 2 3 is smaller than unity 
and therefore the role of matter effects in the P uu oscil- 
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chosen to be 823 = 40°. An overall 10% systematic error has been included in the \ 2 analysis. 
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FIG. 7: Left (right) panel: the different lines depict the 68% and 90% CL (2 d.o.f.) hierarchy resolution using the muon-like 
contained events in a detector with 100 Mt yr exposure, for normal (inverted) hierarchies, and the atmospheric mixing angle 
chosen to be 823 = 50°. An overall 10% systematic error has been included in the \ 2 analysis. 



lation probability is smaller (see Eq. I|15|> ). 

For the inverted hierarchy case, the best sensitivities 
are reached when 023 = 40°. This can be understood in 
terms of the muon (anti) neutrino disappearance prob- 
ability in the presence of matter effects and non 
negligible solar effects. If one performs an expansion 
up to second order in the small parameters Am^/Am^ 
and sin 0i3, there are two terms in the muon (anti) neu- 
trino disappearance probability equation which depend 
on cos 5. One term is proportional to sin 2023 cos <5, the 



other one is proportional to cos 2023 cos 8, and therefore 
it is only different from zero for 023 7^ 45°, changing 
its sign accordingly to the octant in which 023 lies. For 
the muon antineutrino channel and inverted hierarchy, 
this term, i.e. the term proportional to cos 2023 cos 8 is 
the dominant one among the two terms proportional to 
cos 8 in the muon (anti) neutrino disappearance proba- 
bility equation. For 023 < 45°(> 45°), the term propor- 
tional to cos 2023 cos 8 maximizes (minimizes) the impact 
of matter effects in P uu . The former term is also the 
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one responsible for the change of shape in the sensitivity 
curves for 9 23 = 40° and for 9 23 = 50°. For 6 23 = 40° 
the maximum is reached now at 6 = ±180° and not at 
6 = 0°. For the 623 — 50° case, the sensitivity curve 
bends in the opposite direction, due to the opposite sign 

of COS 26*23- 

In Figure[5]we show the results for the case as in Fig. 03 
but here we also marginalize over the value of lAmjjjJ in 
the allowed range. The ability to distinguish the hierar- 
chy appears significantly decreased, after using the same 
2 d.o.f statistics to define CL. It is however important 
to note that we expect the value of jAm^l to be known 
with much higher precision on the time scale relevamt for 
the IceCube deep core experiments. 

VI. DISCUSSION AND CONCLUSIONS 

The IceCube detector and its proposed deep core array 
provide a great opportunity for studies of atmospheric 
neutrinos. Being the largest existing neutrino detector, 
it will accumulate a huge number of atmospheric neutrino 
events over an enormous energy range, thus allowing for 
detailed studies of oscillation physics, Earth density, at- 
mospheric neutrino fluxes and new physics. In order to 
extract all this information it is necessary to use energy 
and angular distribution information, as well as flavour 
composition, all possible to obtain with the IceCube de- 
tector. 

Qualitatively, there are three main energy intervals and 
three main angular regions which are sensitive to different 
types of physics. 

At very high energies, above 10 TeV, neutrino inter- 
action cross-sections become high enough that neutrinos 
going through the Earth start getting attenuated. This 
effect is sensitive to neutrino interaction cross-sections 
and to the density profile of the Earth. A measurement 
of the neutrino flux at these energies can provide a de- 
termination of the Earth density [29| which can be used 
as an experimental input to our analysis, instead of the 
PREM predictions. 

The "intermediate" energy region, between 50 GeV 
and 1 TeV can provide good information about the at- 
mospheric neutrino flux, which can be used to improve 
the uncertainties in the simulated atmospheric neutrino 
fluxes. 

In our paper we concentrated on the "low" energy 
region, below about 40 GeV, where neutrino oscilla- 
tion effects can be significant. Matter effects inside the 
Earth are very important in this energy range and for 
non-zero values of #13 resonance effects can strongly en- 



hance/reduce oscillation probabilities. Straight-up-going 
neutrinos (eg < —0.7) pass through the core of the Earth 
and are most sensitive to resonant matter oscillations and 
thus to all sub-dominant neutrino oscillation effects (#13 , 
mass hierarchy, CP violation). Up-going neutrinos at 
shallower angles are still sensitive to the "main" oscilla- 
tion effects (Am 2 ^, #23), while the other, sub-dominant 
contributions become smaller, due to the lower matter 
densities and shorter pathlengths. 

In the present study we have exploited muon neutrinos 
and antincutrinos with energies in the 10 — 30 GeV range 
which have crossed the Earth (i.e. eg < —0.7). Using the 
/1 like contained events for an exposure of 100 Mton yr 
in the proposed Icecube deep core ice Cherenkov detec- 
tor, the neutrino mass hierarchy could be extracted at 
the 90% CL if sin 2 26>i 3 > 0.02 (sin 2 26 13 > 0.14 when a 
10% systematic error is included in the analysis) regard- 
less of the value of the CP violating phase S. In addi- 
tion, downgoing neutrinos in the deep core will provide 
a measurement of the atmospheric neutrino flux, helping 
enormously in diminishing the systematic uncertainties 

The Icecube deep core array, with muon energy detec- 
tion threshold of ~ 5 GeV could provide the first mea- 
surement of the neutrino mass hierarchy if sin 2 29i 3 is 
close to the present upper limit. This would be possible 
with an exposure of 100 Mt yr and if systematic errors 
below 10% can be achieved. The next generation of long 
baseline experiments would then have to concentrate only 
on the extraction of the CP violating phase S. 
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